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I. INTRODUCTION 



Gravitational theory allows a conserved current J^j^ to be formed from the Einstein 
tensor with a symmetry generated by Killing vector k^: 

ja _ na 1.13 

The current J^j^ is conserved since the Einstein tensor is divergence-free and /c^ satisfies 
Killing's equation. In a similar fashion we define 

"^BACH •= B'^pi^ (1) 

where -8^/3 is the Bach tensor, a symmetric, divergence-free, trace-free tensor [|l|, 0, [Q, 
is a generator of conformal maps which include Killing symmetry and homothety 
as special cases, ^bach is conserved because the Bach tensor is divergence-free and the 
symmetrized covariant derivative of the conformal generator S,^ is ^{a;P) = '^^{x'^)gaf3- Since 
the Bach tensor is trace-free it follows that 3^(3^^"'^^ vanishes. 

The Bach tensor is a second derivative of the Weyl tensor Ca^ipu plus a Ricci tensor term 

Bal3 = '^^'^^Ca^py + -R^^Ca^pu- (2) 

For metrics with non-zero Bach tensor, we can think of Bapig) = K,weyiSai3 as the field 
equations of conformal Weyl gravity 0, with symmetric Sap analogous to Tap of Einstein's 
theory. 

There exist superpotentials for Jgjj^ which lead to an Einstein mass by integrating the 
superpotential over two-surface S^. Three of the superpotentials do not require a coordinate 
choice: the Komar potential 0, the Penrose-Goldberg potential [0], and the Taub potential 
ip. However, each of those three superpotentials is restricted to special systems. For more 
general potentials derived from the Hilbert action and pseudotensors, the Einstein energy is 
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known to be ambiguously defined because it depends on the choice of background coordi- 
nates. Nevertheless, one may use the two-surface construction to define sensible notions of 
total gravitational energy. Indeed, for asymptotic flatness, say, towards future null infinity 
X"*", 5*^ tends to a "nice" cut of X"*", a round two-sphere. The (suitably unique) choice of 
asymptotically Cartesean coordinates ensures that the Einstein energy agrees with the ac- 
cepted Bondi-Sachs notion of total energy ^ within some quasilocal (i.e. finite) two-surface 
S"^. One of the first fundamental expositions of quasilocal mass (and angular momentum) 
was given by Penrose |]lO|. A more recent and well known treatment of quasilocal energy is 



the work of Brown et al 111 



In this work we introduce the conserved Bach current as a tool, in general relativity, to 
study collapsing star models. In particular, we use the properties of the Bach tensor and 
Bach current to study collapsing spherical fluid interiors matched to a Vaidya exterior and 
analyze the power flow. generates a matter symmetry in the Vaidya exterior, weaker than 
conformal symmetry, but C,(a;i3) annihilates the Bach tensor and so the current is conserved. 
The volume integral of Jbacu l^cids to a two-surface integral which provides a quasilocal 
value for the Vaidya power gradient. Covariant thermodynamics of irreversible systems is 
used to interpret the Bach current for the Vaidya exterior. The Vaidya power gradient is 
related to the rate of entropy increase and the number of states the energy resides in. 

We can view a star model as a single state drawn from a microcanonical ensemble. The 
initial values of all the observables are at their equilibrium values (say, all zero), then, as a 
function of time, they will not remain exactly zero but will execute small fluctuations about 
their zero equilibrium values. The Vaidya exterior dissipates energy by sending "short- 
wavelength photons" to null infinity. It is the Bach current and the rate of entropy increase 
that allows one, via the fluctuation-dissipation theorem, to measure the rate at which a 
deviation from equilibrium decays. 
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CONVENTIONS AND NOTATION 

In this work Greek indices range over (0,1,2,3) = {u,r,i),ip). Sign conventions are 
2^iy;[a/3] = ^^(3 = R^ajSy Overdots abbreviate d/du and primes abbreviate 

d/dr. The metric signature is (+,-,-,-) and the field equations are G^y = —SnT^iy. 

II. SOME PROPERTIES OF THE BACH TENSOR 



1. A variational derivative of the Lagrangean density C = ^f—gCajs-ysC"'^'^^ yields the Bach 
tensor 

and the Lagrangean provides the field equation 

Bapig) = 0. (4) 

Solutions of Eq.@ are vacuum solutions of conformal Weyl gravity P). 

2. Kozameh et al [0 proved that any solution of Eq.(^) which is also conformal to C spaces 
is conformal to an Einstein metric (Einstein metrics have Ricci tensor proportional to the 
metric). For an example of a metric which satisfies Bajsig) = but is not an Einstein metric 
see Nurowski and Plebanski p]2| . 

3. Vacuum spacetimes, conformally fiat spacetimes, and Einstein spacetimes have vanishing 
Bach tensor. 

4. For Pfj^i, = R^j_y — {R/6)gfj,iy, the Bianchi identities allow the alternate expression 

Bo^p = V^Vio.P,]p + \p^'C^,p,. (5) 

5. Under conformal map cjap = VL^gap, the Bach tensor transforms as Bap = Bap. 

HI. THE VAIDYA SPACETIME 
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A. Vaidya Interiors 



The first Vaidya interior was constructed by Synge In that work an incoherent 

shell of radiation is joined to Vaidya. Santos \T^\ joined a collapsing fluid with radial heat 
flow to a Vaidya exterior. Glass [|T3| studied a two-fluid interior, a collapsing perfect fluid 
with an outgoing null fluid matched to a Vaidya exterior. Continuity of the first and second 
fundamental forms expressed conservation of momentum across the boundary and revealed 
the boundary to be the perfect fluid zero-pressure surface. This is the interior we have in 
mind to model the energy flow considered here. As the perfect fluid collapses, energy is 
removed by the null fluid. In the Vaidya exterior, the outgoing short-wavelength photons 
carry energy to future null inflnity. 



B. Vaidya Exterior 

The Vaidya metric in outgoing null coordinates, with arbitrary mass function m{u) and 



A := 1 — 2m{u)/r, is given by 

glfdx'^dx'^ = Adu^ + 2dudr - r'^{dd^ + sin^ d d^^). (6) 

We use a Newman-Penrose null tetrad: 

l°'da = dr, ladx^ = du, (7a) 

n"a„ = du- {A/2)dr, n^dx" = dr + {A/2)du, (7b) 

m'^da = {r-^/V2){d^ + ^9^)^ niadx'' = -{r/V2){d^ + isin^ dip), (7c) 

with non-zero spin coefficients 

p = -1/r, /i = -A/(2r), (8a) 

a = -coW/{2V2r) = -p, (8b) 

7 = A' /A. (8c) 

The Weyl tensor is given in terms of a basis set of anti self-dual bivectors (^1^2 = —m{u)/r^) 
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= [m{u)/r^]{U'"'V^'' + V'^'U^'' + M'^^M''^ + c.c), (9) 

where 

^ 2m["n^l, M"^^ = - 2m^'^Th^\ = 2l^''m^l (10) 

The Vaidya spacetime is Petrov type D and has the invariant relation 

RaP^.R''^'"' = C^P^.C^^^^" = 48V1/2. (11) 

The Ricci and energy-momentum tensors are 

R^, = -SvrT^, = ^—j l^l,. (12) 

The standard interpretation of the Vaidya energy flow was given by Lindquist et al |T6[ with 
respect to an observer at rest at future null infinity whose worldline tangent is 

u^d, = A-^l^du = A-^l\it + \AF)d^. (13) 

The local energy flux / is 

^^^^"^^ = - (t^) 1 o\ ^, ^ (14) 

and the total luminosity seen at infinity is L^{u) = j!™^ inr'^f = —m. 

Substitution of the Weyl tensor Eq.(^ and the Ricci tensor Eq. (p!2D into Eq.(|^) yields 
the Bach tensor 

.2m,, , .2m, 



Bap{g^'"^) = {^)Lh + {^){Lnp + njp + m^mp + m^nip). (15) 



The independent eigenvectors of Bap{g^^'^) are dr = l°'da, = {r/\/2){nf + 171°") da, and 
= —i{r / ^/2)smi9{m°' — fff)da- 

IV. VAIDYA CONSERVED CURRENT 

The generator of time translations is = ^"<9q,. The time dependence of m{u) prevents 
from being a Killing symmetry but generates a matter symmetry since C^g^l^'^ = 
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—rRa(3- The Bach current Jbach — -^"^^/S metric g^"''^ will be labelled Jyad- The four- 
current is conserved since Bap^^"'^^ = B^fsl^l^ = 0. Stronger symmetry can be imposed by 
restricting the mass function. With a linear mass function, the Vaidya metric has been used 
to study self-similar collapse il^. This restriction was anticipated by the author |TH] for 



m(ui 



tuq + rriiu and homothetic symmetry ^"9^ = (mq + u)du + Or where jC^g^^ = '^Uc$ ■ 



With = (yl/2)/" + n"", the conserved four-current is 

r.,. - (^) + + A%) r. (16) 

The four-current is more transparent when written in terms of curve tangents 

•'?.A=(^)a.+ (^)a, (17) 

and is similar to a conserved Maxwell current \/—gJMax — y/~9{Pi J)- By analogy we 
call 2m/r^ the Bach "charge density". Examining ^/—gJvad ^^^^ one that 2m/r^ is the 
flow across u = const null three-surfaces and 2m /r is the flow across r = const timelike 
three-surfaces. This would be equally apparent upon integrating da{\/—gJvad) — '^^^^ ^ 
four-volume bounded by four such surfaces. 

A potential exists for Jy^^, namely the bivector density 



Uf., = v^(9)4/["n^l (18) 



with \/—gJvad = ^iJ^vad- We integrate over a. u = const null three-surface J\f: 

J V^JvaddSa = / Uf^,dS^p (19) 
AT dAf 

where dSa = ladrd-ddip and dM is a closed two-surface bounding M with dS^p = ([anpid-ddip. 
The integral of the Bach "charge density" over dAf provides the quasilocal Bach charge 
within dJV, which here is the Vaidya luminosity divided by luminosity distance (and so we 
call it "luminosity gradient") 

/(-)sin^?4/['^n^]rf5,« = -Stt-, (20) 
J r r 

where r ranges from the spacelike surface r = 2m to future null infinity. 



V. THERMODYNAMIC INTERPRETATION 



The Eckart theory of dissipative fluids is the simplest relativistic generalization of Navier- 
Stokes theory, but Eckart's theory suffers from a lack of stable solutions and acausal propa- 
gation of perturbations. We will follow Calzetta |TU| and use the framework of more causal 



and stable theories such as the Israel-Stewart 'second order' type theory or the Geroch- 
Lindblom 'divergence type' description of a relativistic real fluid. A covariant theory has 
the following rules: 

(a) Intensive quantities {T,p,fi) are associated with scalars, which represent the value of 
the quantity at a given event, as measured by an observer at rest with respect to the fluid. 

(b) Conjugate extensive quantities {S, V, N) are associated with vector currents 
{S°',u°', N°'). If any of the extensive currents X° is conserved, then = 0. The quantity 
■u" associated with the volume is the fluid unit four-velocity. 

(c) Energy and momentum are combined into a single extensive quantity and associated 
with the tensor T°'^. The energy current is = T'^^up. 

The entropy current 5*" is given by 

which is rewritten as 

5° = T"^/5^ + $° - aAT" (22) 

with affinity a = fi/T, thermodynamic potential $° = pf3°', and inverse temperature vector 
pa _ {^ot^rp^ jya _ where n is the particle number density seen by a comoving observer. 

Suppose the fluid departs from equilibrium by a fluctuation 6N°', 6T°'^, consistent with 
the conservation laws but otherwise arbitrary. Then the change in entropy production is 

5^", = -«,„5iV" + /5ft,5T"^. (23) 



For a true equilibrium state the entropy must be stationary ||2^, and so a^a = f3{i3;a) = 
Thus the affinity must be constant, and the inverse temperature vector must be Killing. 
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If generates a timelike conformal symmetry then one has thermal equilibrium with 
inverse temperature jS'^ = ^^/T as a conformal Killing vector for trace-free 5T"^. In the 
exterior Vaidya region we have a weaker matter symmetry C^g^^'^ = —rRap and so SS'^^ 7^ 0. 
The physical dimensions of S are energy/ (sec °K) and, reasoning dimensionally, 



S (24) 



in units with ks = G = c = 1. If one examines the terms of Eq.(^) its clear that and 
A^" lie along u^da = A~^^'^du- T°'^/3f^ has components along and orthogonal to u". We can 
thus write 

where P°'Ua = 0. The divergence of S°' is 

SI = SA-'/^[S/S + A-\m/r)] + (25) 

wherein the Vaidya luminosity gradient appears. The logarithm of S measures the number 
of available states and we see here that rh/r provides a measure of the rate of change in the 
number of states. 



VI. DISCUSSION 

Luminosity is one of the key observations made of astrophysical objects (a mass- 
luminosity relation is an important feature of main sequence stars). Use of the Bach current 
to compute luminosity gradients provides a covariant method for obtaining an astrophysi- 
cally important quantity. We have identified the luminosity gradient with entropy production 
in the Vaidya exterior. The Vaidya exterior dissipates energy by sending " short- wavelength 
photons" to null infinity. For an equilibrium ensemble of star models, it is the Bach current 
and the rate of entropy increase that allows use of the fluctuation-dissipation theorem to 
measure the rate at which an initial deviation from equilibrium decays 
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A full development of dissipative thermodynamics using the Bach current requires ana- 
lyzing the collapsing Vaidya interior, which we leave for future work. The Vaidya spacetime 
is only a first approximation to interesting systems. 
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